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Chapter 11 Integration on manifolds 


We are now almost ready for our concluding chapter on the great theorems of classical 
vector calculus, the theorems of Green and Gauss and Stokes. 

The final thing we need to understand is the correct procedure for integrating over a 
manifold. Of course, manifolds are typically curved objects, so there are significant issues 
here that we did not have to face in dealing with integration over (flat) Euclidean space. Even 
though we have been able to use integration to compute such things as the volume of a ball 
in R”, yet the integration involved is essentially flat. This is reflected in a formula such as 





dvol, = dx, dt2...dXp. 











As a specific example, we have for the ball B(0,r) in R®, 





4 
vol3(B(0,r)) = rue 
but we have not yet even defined what we should mean by the area of the sphere, 
volo(S(0,7)). 


It will in fact turn out that there is a very reasonable way to accomplish this task, and it will 
be based on our knowledge of flat integration. 

We choose to employ the parametric presentation of manifolds as discussed thoroughly 
in Section 6E. It is in this context that we shall see how to define volume and integration. 
Though our definition is to be based on a particular choice of parameters, we shall be able to 
prove that the integration we define is actually invariant under a change of parameters and is 
thus intrinsic to the manifold. 

We denote by M an m-dimensional manifold contained in R”. In what follows we shall 
first investigate how to define the m-dimensional volume of subsets of MW, and after that we 
shall easily define integrals of real-valued functions defined on M. 

We begin with sort of a warm up case, that of one-dimensional manifolds in R”. 

















A. The one-dimensional case 

We assume in this section that m = 1, so we are dealing essentially with a curve M C R”. 
We effectively already know exactly how to deal with this case, thanks to Section 2B. 

If M is represented parametrically as the image of a one-to-one function of class Ct, 





(a,b) > M, 
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then the length of M/ is given as 


vol|(M) = / \| F(t) ||at. 





Of course, M Cc R” and F"(t) € R” and ||F"(t)|| is the norm of the vector F’(t). 














We really need to say no more about this definition except to remark that if we imagine the 
interval (a,b) to be partitioned into small pieces (a, 3), then F’ sends (a, 3) to an approximate 
interval in R” whose length should be approximately || F’(a)||(@G—a@). Thus || F’(a)|] represents 
a scale factor relating the parameter length G—a to a length in R”. A useful way to represent 
this definition is 


























dvol, = ||F’(t)||dt. 


As we discussed in Section 2B, the definition of vol;(/) is independent of the choice of 
parametric representation of M/. 











PROBLEM 11-1. Suppose M C R? is a curve described in polar coordinates by an 
equation r = g(0), where a < @ < b. Show that the length of M is 


/ JaO? + g(O)2d0. 


(This formula is usually written as 


dr\? 
iS ieee 2 
/ (F) +r? dd 


r J (dr)? + UBF 





or even as 
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PROBLEM 11-2. A certain curve in R? is described in polar coordinates by the 


equation 
= */ cos 20. 


This “figure eight” curve is called a lemniscate. (See Problem 2-6.) Show that its length 
is equal to 


ra) _ rt) 
VETS) ~ an 





r 


PROBLEM 11-3. Show that 





n/2 d 
i Vi+sin20 


(HINT: substitute « = sin@ and then x = t'/*.) 


PROBLEM 11-4. Show that 
7/2 sin? g 


ff = 
0 V1+sin?6 


PROBLEM 11-5. Combine the preceding two problems to show that 


[vir svt = Vt P (3) gee: 


a NG) P (3) 


PROBLEM 11-6. The two cylinders x? + z? = 1 and x?+ y? = 1 intersect in a certain 
curve. Find its length. 


(Answer: 8 times the answer to the preceding problem.) 





Finally, we easily extend the definition of the length of WV to the more general definition 
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of integration of a function over WM. Suppose then that M _, R. Then we write 





b 
[ favor = f POF ole 
M a 
In this particular case this is often called an arc-length integral, and is written 


| fds or fdé. 
M M 


A frequent application of these integrals is the calculation of the average value of f over 
M, defined to be 


= 1 
= Se dvol,. 
i vol, (MM) we ved 
Notice the easy facts 
ft+g = f+5, 
cf =cf_ ifcis a constant, 
1 ee 











DEFINITION. The centroid of the curve M is the point  € R” whose coordinates are 
given by 





t;, = average of x; over M 


1 
— ,dvol,. 
vol, (M) [ sadvot 


In other words, z is the average position in R” of points in M. Note that it is not to be 
expected that x necessarily belongs to M. 














PROBLEM 11-7. Calculate the centroid of a semicircle. 


PROBLEM 11-8. Calculate the centroid of the “right-half” of the lemniscate of 
Problem 11-2. 


(Answer: #=4,./n/T? (1/4), y=0.) 


B. The general case 
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Amazingly, what we have accomplished in the one-dimensional case generalizes almost 
immediately to m dimensions. We thus suppose WM C R” is an m-dimensional manifold. We 
also suppose that M is given by a parametric presentation of class C', 














A->R", F(A)=M 








where A is a contented subset of R”. For short we may write x = F(t), so that the coordinates 
of x are given real-valued functions of t = (t),...,tm). In our discussion we must have 
1<m <n; the case m = n is tantamount to a discussion of integration over R” itself. 

In every situation our goal is the following: to use the calculus associated with F' to find 
a scale factor Z(t) which converts the m-dimensional “infinitesimal” volume dt = dt, --+ dtm 
of the parameter space to the corresponding m-dimensional volume of /. We shall write this 
correspondence in the form 




















avol«, = J (dt. 








We have chosen the notation 7(t) to remind ourselves that this scale factor is like a Jacobian 
determinant. In this and other formulas of a similar nature we shall be very insistent upon 
gaining an understanding of the Jacobian Z(t). 

As in the one-dimensional case, what we really need to do is see how F' transforms small 
rectangles in A and in particular what scale factor to multiply by to compute the approximate 
m-dimensional volume of the image of the rectangle. As usual, let €1,...,€,, denote the unit 
coordinate vectors in R™. Then consider a special rectangle 








nN 
t+e,e, 


[é1,t1 + €1] <--- X [tm tm + Em, 





Nn 
t t+ €,€, 


where ¢ is a vertex and the edge lengths are €1,...,€m. We think of the €;’s as small positive 
numbers. 

The corresponding image of this rectangle with respect to F is approximately an m- 
dimensional parallelogram with one vertex F(t) and edges given by 


OF 
F(t) Ea 
1 


Bp, Or ema): Ft) + a5, 


F(t) 


6 Chapter 11 


Of course this flat parallelogram does not lie in the manifold M. However, it is tangent to 
at the point F'(t). 


R” R” 











ty 


We know how to compute the m-dimensional volume of this parallelogram: from Section 8C 
it is 








/Gram(€,0F/0t1,...,€mOF/Otm) = €1...€m\/Gram(OF/Ot,...,OF /Otm). 


Since the factor €,...€;, is the volume of the rectangle in the parameter space, we see imme- 
diately that the scale factor for computing the volume of / is precisely 





/Gram(OF/0t,,...,0F/Otm). 


Since we think of e; as dt; in some sense, and thus €, --- €,,, is like dt, ---dt,, = dt, we are thus 
led to the 


DEFINITION. Given the C! parametric presentation of the m-dimensional manifold M Cc 
R” 





A—>M, 











where A is a contented subset of R”, 








voln(M) = | /Gram(OF/dt,,...,OF /Otm)dt; ... dtm. 
A 


More generally, if f is a real-valued function defined on M, its integral over M is given by 





| fav0ln = [ 1) VGram(OF[80,... OF [Atha 


Integration on manifolds 7 


We summarize both these equations by simply writing 





dvolm = /Gram(OF/0t1,...,OF /Otm) dt. 


This equation gives the scale factor 7(t) explicitly in this general setting. 
There is a nice alternate description of the scale factor described above. Namely, use the 
Jacobian matrix 





ar aFy 
Oty “'° Otm 

DF =| : : 
OF n OF n 
Oty **"°  Otm 


(see Section 2H). This Jacobian matrix is an n X m matrix, and the matrix product 
(DF) DF 


is precisely the m x m Gram matrix (OF /0t; e OF /Ot;). Thus the scale factor in the above 
definition can also be written in the form 


I(t) = /det(DF)'DF. 


EXAMPLES. For the case of curves, m = 1, this is the same definition we used before, as 
Gram(OF'/0t) = Gram(F"(t)) = F’(t) e F’(t), so the scale factor is || F'’(t)]]. 
For the case of 2-dimensional manifolds, we have 





voly(M) = | JOF/ati 2 ||OF /ate||? — (OF dt, e OF /ate)2dt. 











Let us examine the sphere $(0,a) of radius a in R?. Let us use spherical coordinates in the 
usual way, so that 





F(y,@) = (asinycos@, asinysin@, acosy). 


Then 


OF /Op a(cosycos@, cosysin#é, —siny), 
OF /00 = asiny(—sin@, cosé, 0). 


The scale factor is therefore 





Ja2(asiny)? —0 = a’sing. 
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Thus for this sphere we have 


dvoly = a’sin pddd. 


In particular, the area of the sphere is 


20 wT 
volo(S(0,a)) = | | a’ sin pdydd 
0 Jo 


= Ana’, 


a familiar result to be sure. 





The example of a right circular cylinder in R® is even easier, but it is so important we 
should pause to consider it. Let us use the cylinder of radius a with axis given as the z-axis: 


It is convenient to use cylindrical coordinates (no surprise!), which we call 0, t: 


x =acosé, 
y =asind, 
2 =: 


We thus have 


F(6,t) = (acos@, asin 6, t); 


OF : : 
a a(— sin 6, cos 6,0); 
OF 
OL —= (0,0;.1); 

SL) Se 


dvoly = adédt. 
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PROBLEM 11-9. Consider the right circular cone as in Problem 6-29: 


h 
z 


z= cota,/r?+y?. 





In terms of cylindrical coordinates 6, z show that for this manifold 








dvolg = ztanasecadédz. 








Then show that the area of a right circular cone whose base has circumference c and 
whose slant height is 0 equals scl. 


Just as we have done in the one-dimensional case of the preceding section, we need to check 
that our definition is independent of the choice of parametrization of M. In other words, we 
must verify that if we change parameters we do not change the value of the integral. To be 
sure, our geometric intuition tells us that all is well, but we should do the work to certify the 
formulas. 


Thus we suppose that M is parametrized in two manners: 
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D> .,. 


. m : n 
coordinates ¢ R' coordinates s R 














There is a parameter change involved here, namely 
h(s) = F-'(G(s)). 
That is, h = F~'o G; we may represent 
i = fils), 
and h is a C' function from B onto A and has a C’ inverse. 


THEOREM. In the above situation the integral 


fdvoly 


M 


has the same value whether computed using F or G. 


PROOF. By now we realize that this proof will probably amount to no more than an 
automatic calculation. And so it does. The scale factor in using G is the square root of 
det(DG)'DG = det D(F oh)'D(F oh) 
= det (((DF)oh)Dh)' ((DF)0h)Dh_ (chain rule) 
= det(Dh)' ((DF) oh) ((DF)oh)Dh 
det(Dh)! det[((DF) 0 h)' ((DF) 0 h)|det Dh 
(det(Dh))*|[det(DF)'DF] 0 h. 
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Thus we obtain 





[ f(G(s)),/det(DG(s)'DG(s)ds 


= [ Fes) lace Dis) det(DF)'DF(h(s))ds 





Z | f(F(t))\/det DF) DF (hat. 
A 


This is the desired equation. Notice that the last step used the change of variables formula 
for integrals over R™, t = h(s) and dt = | det h’(s)|ds, which we introduced in Section 10G. 

QED 

In summary, it should be stressed that the thing which shows that our definition of inte- 

gration over manifolds makes sense is precisely the change of variables formula of Section 10G. 














PROBLEM 11-10. The flat torus in R* is a nice example of a two-dimensional 
manifold. See Example 3 of Section 5G for the definition. Calculate its two-dimensional 
volume (its area). 





C. Hypermanifolds 
We now examine the important special case of a manifold M C R” of dimension n—1. The 
scale factor for volume calculations has a particularly nice form. We first give a calculation: 





LEMMA. det (6; + aaj) =1+ D5 aj. 
j 


PROOF. This is a special case of Problem 3-41; nevertheless we provide a proof here. 
Suppose we are dealing with an n x n matrix. Then its j column is €; +a,;a, where a is the 
column vector with coordinates a,,...,@n. The linearity property of det (Section 3F), implies 
that the determinant we want, 


det (é, + a1aQ €2 +Q9Q@ °*:: En + AQ) ; 
is equal to the sum of 2” terms, each term being the determinant of a matrix with either 


é; or aja in the j™ column. If a scalar multiple of a appears in at least two columns, the 
corresponding determinant is zero. Thus only n of the terms of this sum survive, and we 


12 
obtain for the above determinant 


det I+ S © det (é: “++ Gj +++ En) 
j=l 

2a eae) 
j=l 


= 1+ 5 a;det (é ers aj€; ns Ey) 
j=l 


= |+ Sr 
j=l 


PROBLEM 11-11. For the vectors €; + a,;é,, 1 < 7 <n -—1, show that 


n—-1 
P , ‘ 2 
Gram(é, + a1€n,.--,€n-1 + @n—1€n) = 1+ y Qj. 
j=l 
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QED 


Now suppose that the hypermanifold IM Cc R” is presented as the graph of a function 




















A —+ R, where A C R"~!. Then we obtain an explicit parametrization of M with A = R”, 








where 
FB Gigs ney Qt) = ign nng Pacts Py <n <5 ncd)): 
Therefore ae 5 
mn Y. ‘ 
z~— = + 7—€n;, lagen 1, 
Ox; "3 Osc; J a 


Problem 11-11 immediately gives 


n—-1 
Gram(OF/0x1,...,0F/O%,-1) = 1+ 5 > (Op/dx;)’. 
j=l 


Thus we can compute volume integrals over MV in terms of the coordinates 7}, . 





n-1 
dvol,_; = ,/1+ (Op) Oms)* diy xn Chn4: 


j=l 





++5Un-1; by 
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PROBLEM 11-12. Compute the area of the portion of the surface z = xy contained 
inside the cylinder x? + y? < 3. 


(Answer: 147/3.) 
PROBLEM 11-13. Compute the area of the portion of the paraboloid 2z = x? + y? 
contained inside the cylinder x? + y? < 3. 


We can now immediately generalize this result to hypermanifolds presented implicity. 


THEOREM. Suppose M C R” is a hypermanifold presented as the level set 





n 





where g is a function on R” of class C' satisfying Vg #0 on M. Suppose for instance that 
Og/0x, #0. Then integration over M can be calculated by means of the formula 


I|Vall 
dvol,1. = —— dz, ... d&y_1. 
|0g/Ox| 
PROOF. The implicit function theorem guarantees the local existence of a function R"~! > 








R such that 





g(x) 0 > tp O Hi peep ind). 





Furthermore, differentiating the equation 


Ol Ws oss 5 Ol Pig en ss e)) S 
gives the formula 
0g/0x; 
0g/02n 


Thus the scale factor found above in terms of y can be manipulated to give 





0/02; = — 








A 


n— 


(Op/Ox;)? 


1 


1+ 09/0%,)? | (Og 0%)" 


J 





1 


n—-1 
= ——— ||(0g/0z,)?+ >) (0g/0z;) 
gO] | (29/Fr0)? + 2 (O9/2e:) 


14 


EXAMPLE. Let the manifold be the sphere S(0,a) C 
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QED 











R". This is an (n — 1)-dimensional 





manifold and is described implicitly by the equation 


2 


0. 


Iz? —@ 


For the function we have 


Vg 


2%. 09/08, = 2a: 


In order to cope with 0g/0x,, 4 0 we shall find the volume of the upper hemisphere x, > 0 


and then double the answer. Thus 





2 
voln_1(S(0,a)) = 2. f Gal dx, ... dtp 
=. 2 = dx, Pere ALy—1 





dx, dais Onc 





In 
D / 7 
Ja? —_— a — ee ee: 


where the integration region is the ball 77 +---+22_, <a? in 


therefore obtain 


vol,_1(S(0, a)) 2( 


Substitute r = av/t to obtain 








" a 
n—1)o,_, | ———= 
if Ver 








2 


ol 


1 





R”-!. From Problem 10-49 we 











r 2dr. 





1 
1 
vol, S(O) = Mons | A heat hat 
5 vio 

(n—1) aA Lge 

= (n—1)a,_14 
: 5 LSE 
n—-l1 1 
— (n= TaysaB ( 5 5) 
2h 
ee) l (3) 
Dqpr/2 fe 
lr (5) 
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We record this result for later reference: 





Qnr/? 
vol,—1(.9(0,1)) = T(n/2) = Yi, 














For instance, 


vol2(S(0,1)) = 47, 
vol(S(0,1)) = 2n?, 


vol4(S(0, 1)) 





EXAMPLE. Here we simply reconsider the sphere $(0, a) in R°. In Section B we calculated 
the area formula 
dvolz = a’ sin pdyd 


in terms of the usual spherical coordinates. It is quite interesting to think of this result also 
in terms of cylindrical coordinates for the sphere: x? + y? = a? — z?, so that 


x = Va? — 2? cos8, 
y = vVa?— 2?sind, 


z = &Z. 


Of course we assume here that —a < z < a. Then we can calculate 7(z,@) from scratch, or 
we can instead observe that in terms of the angle y we have z = acosy; thus dz = asin pd 
(we’ve included the absolute value). The result is that 


dvoly = adzdé. 


This is a truly fascinating result. The expression for this area is exactly that of the circum- 
scribed cylinder 
ety = al 
An interesting application of this result has to do with constructing maps of the sphere 
onto planes. If we project the sphere onto the cylinder by moving in horizontal straight lines 
from the z-axis, then area is preserved! The precise formula for this map is 


F( ; ( ax ay ) 
ZY, 2) = $= 8 ee. 
2 Ve@-2’ J#@—2 
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Its inverse is 





x 
F(¢,9,2) = (vo? = =, vo? = 272, 2). 
a a 


this is an amoeba on the sphere 


eS this is its projection onto a cylinder 











Then the map onto a plane is obtained by cutting the cylinder along a generator and unrolling 
it. 

An interesting and curious consequence comes from the above observation. Namely, a zone 
on the sphere is defined as a set of the form 


{(z,y,2) | a +y +2? = ar, 2 2s Bo}: 


Here —a < 2 < 2 <a. The area of this zone is precisely 27a(z2. — 2), as that is the 
corresponding area of the circumscribed cylinder. Thus the area of a zone depends only on 
its “height” z 2 — z,;, and not on its location. 


EQUAL wee 


_— AREAS 
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PROBLEM 11-14. Let H be the unit hemisphere in R® described as x? + y? + 2? = 1, 
O<z<1. Let C be the boundary equator «7+ y? =1, z =0. Let 0 <a < % and let M 
be the portion of H which lies over the triangle shown in the figure. That is, in standard 


spherical coordinates M is described by y 
A 


0 
= 


sin y cos 6 





Show that the area of MW equals 





volo(M) = 2a+acosa—T. 


(HINT: do no integration.) 


PROBLEM 11-15. (This problem is based on Problem B3 from the William Lowell 
Putnam Mathematical Competition for 1998.) Continue with the hemisphere H of the 
preceding problem. Let P be a regular n-gon inscribed in the equator C’. Determine the 
surface area of that portion of H lying over the planar region inside P. 


(Answer: 27—n2+n7 cos 7/n) 


PROBLEM 11-16. Find the area of the portion of the surface of the cylinder 2?+y? = 
1 in R® which is contained inside the cylinder x? + z? < 1. 


(Answer: 8) 


PROBLEM 11-17. Find the area of the portion of the surface of the cylinder x? + y? 
1 in R? which is contained in the intersection of the two cylinders x? + z? < 1, y?+272 < 


(Answer: 16—8V2) 


Ue 
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PROBLEM 11-18*. A surface M C R? is described implicitly by the equation 


(t+ yt2) = oy. 


Prove that this is actually a two-dimensional manifold at each point except the origin. 
Compute its area. 


(Answer: 12/2) 


PROBLEM 11-19. Consider these two surfaces in R*: the unit sphere x? +y?+2? =1 
and the right circular cylinder x? + y? = y. 


a. Find the area of the portion of the unit sphere that is inside the cylinder. 
b. Find the area of the portion of the cylinder that is inside the sphere. 


PROBLEM 11-20. Let M be the “elliptical cone” in R® defined by the equation 





and the inequality 0 < z < h. Show that the centroid of M is (0,0,2h/3). (You will not 
be able to calculate the area of M.) 





D. The Cauchy-Binet determinant theorem 
This section is not very important for our development of integration, but is so fascinating 
that it should be included. 
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We start with an interpretation of the formula from Section 7B: for any two vectors a, 
be R" 














Gram(a,b) = S (aid; — a;b;)?. 
i<j 
We know of course that the left side equals the square of the area of the parallelogram P with 
vertices 0, a, b, and a+b (Section 8C). Now consider the orthogonal projection of P onto the 
coordinate space R*, namely the x; — x; plane. This projection P(i,7) is the parallelogram 
with vertices 





(0,0), (@4,:05), (6y;.05), and (a; + b;, aj + b;). 











The square of its area is of course the corresponding Gram determinant, 
a 
a;+ae ab; + a,b; 2 
a J oe Jed _ b. —q.b: 
det &: taj; BLE ) (debe, 0,)" 
Thus we have the interesting fact that the square of the area of P is equal to the sum of 
the squares of the areas of all the projections P(i,7), for 1 <i <j <n. This is sort of a 
Pythagorean theorem for two-dimensional parallelograms in R”. 

The above result is a special case of the theorem of this section. Before discussing it let us 
see what our special case implies about the change of variables formula for two-dimensional 
manifolds. We use the parametric presentation of MC R”, but instead of naming the presen- 
tation function we simply denote its components as %1,...,%, regarded as functions of ¢1, to. 
Then we have from the above equation 

Ox;/Ot, Ox;,/Ot.\\? 
t = a eee ‘ 
(F( - t2)) 2 (<c G6 Ox; / Oty 


i<j 








In terms of the notation mentioned in Section 10I, this says 


(T(t, t2))? = ye (Beem) 


i<j 


20 Chapter 11 


This makes the volume formula look like a rather natural generalization of the formula for 
changing variables in an ordinary integral over R?: 


al 2 1/2 
dvols = ( A(t, ty) ) | dt,dto. 


i<j 





The above result does not appear to be particularly useful for calculations, but its elegance 
is very interesting. The algebraic generalization we now present is likewise not something we 
shall be using, but we present it here to display a beautiful piece of algebra. 

Instead of m = 2, we now consider the general case of 1 < m < n. For the present 
discussion we shall use the notation u to stand for a generic m-tuple of increasing integers in 
the range [1, n}: 

= (igen g tin). Lh =< ie oo Sy, ~ 
Given an m x n matrix A and an index u, we can construct a (square) m x m matrix A(w) 
by using only the columns of A corresponding to u: 


Qhu, shies Alum 
A(u) = . . 


Amu, +++ Alum 


Likewise, given an n X m matrix B we construct B(u) by selecting only the rows of B corre- 

sponding to u: 

Did! ose. Diam 

Bu) = : : 
ied «ae. Vat 


Here is the amazing result: 


THEOREM. /f1<m<n and A is anm x n matrix and B is ann x m matriz, then 


det AB = S- det A(u) det B(w). 


PROOF. The ij entry of AB is of course 


y Aikdk;; 
k 


Integration on manifolds 21 


where the sum extends for 1 < k < n. Therefore we can write the j** column of AB in the 
form 


(AB) cotumn go » Dis (A eeician k- 
k 


Now we use the linearity of det with respect to its columns, the multilinear property of 
Section 3F. For each column of AB we need a distinct summation index; so in the above 
formula for the j** column we change k to k;. Then we have 


det AB = »; bk1 > + Dkmm det ((A)cot ky °°" (A) col km) « 
i 


In this huge sum each k; runs from 1 to n. 

Now think about the determinants in the sum. If two of the k,’s are equal, the corre- 
sponding matrix has two identical columns so the determinant is zero. Thus we may restrict 
the sum to distinct k,,..., km. These distinct choices, if reordered, correspond to the various 
increasing sets of indices u. For a given u, there are m! distinct orderings k,,...,k,, and each 
determinant in the sum will be +det A(u). Thus we can rewrite the above formula as 


det AB = S~ f(B,u) det A(u), 


where f(B,u) just stands for some number that depends only on B and u (and not on A). 
Now it’s easy to compute each coefficient by just picking a particular A. Given a particular 
u choose A to be that matrix defined by 


Qiu; Oi; 
Oe = 0 TE BA tijxcag tins 
That is, the use column of A is the coordinate vector e; € R", and all other columns of A are 
0. Thus all det A(u’) except u’ = u are zero, and 
det AB = f(B,u) det A(u) 
f(B,u) det I 
On the other hand the 77 entry of AB is 


y Qindej = Gin buggy + +--+ Cin Oumg 
k 





= bug +++ + OimBuns 
=) 


uid 
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That is, AB = B(u). This proves that f(B,u) = det B(u). 


REMARKS. 1. This result in case m = n just asserts that for square matrices det AB = 
det Adet B. Thus we have given a third proof of this crucial property of determinants! 


2. This proof should remind you of the proof of the lemma of Section C. In fact, in that 


notation 
(day + aja;) = B'B, 
where 
1 0 
0 1 
B= eS 
0 O 
ay ag 
is an (n+ 1) x n matrix. Note that 
2] (Clee) 2s i 
1 
B((1,...4-Le+1,....24+1)) = : 
ay 


sO 


det BA vieuyn)) = 1, 


a2 


An, 


ha row omitted 


det B((1,...,.k-1,k+1,...,2+1)) = +ay. 


Thus the lemma is a special case of the theorem. 


3. In general, in case A = B*, we have for m x n matrices with 1 <m <n, 


det AAt = > (det A(u)) 


U 





4. Consider vectors x, y, u, v in R® and define 


He “1 XQ 
YU Ye 
Uz Ui 


U2 V2 
U3 U3 


B 


X3 
¥3 


7 


p) 
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ae tees a . 
yeu yeu 
so the theorem gives 


aor (22% £28) = cee (2%) ae (#2) 
yeu yeu Y Ye U2 V2 

+ det & *) det (‘: a + det (* . det (‘ : 

Yi YB U3 U3 Y2 YB U3 U3 


This is precisely the result of Problem 7-3, the generalization of Lagrange’s identity. 


Then 


E. Miscellaneous applications 


1. Surfaces of revolution 





We are going to examine a special class of surfaces in R® obtained by revolving about an 
axis. We may as well use the z-axis as the axis of revolution. We then start with a curve y of 
class C' lying in the x — z plane, where all points of y have positive x-coordinates. 








We revolve the points of y about the z-axis through a full angle of 360° and thus obtain a 
2-dimensional manifold (a surface) M. We want to investigate integration on M. 

Let y = y(t) be a parametrization, where a < t < b. We assume 7¥ does not intersect itself, 
except that perhaps y is a closed curve (y(a) = 7(b)). We shall designate the coordinates of 
7 as 


Thus, x(t) > 0. 
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Let 6 denote the angle of rotation about the z-axis. Then M can be parametrized as 
follows: 
F(t,0) = (a(t)cos@, x(t)sin@, 2(t)). 


This is a representation using cylindrical coordinates, where the distance from the axis is x(t). 
Then we compute easily 


OF = (x'(H)cos8, a'(t)sind, 2'(0), 
_ = x(t)(—sin6, cos6, 0). 


These tangent vectors are orthogonal, so the corresponding Gram determinant is 


a 


Thus the Jacobian factor for dvolz is 


OF 


06 




















= (a'(t)? + 2(t)*) a(t)’. 


FEC) = 2s ey + et. 


Hence we obtain for instance 


voly(M 


[ a (t) Ja"? + 2 (pdt 
Qn i a(t) /a! (t)? + 2/(t)2dt. 


area(M) 


Another way of saying this is 


area(M) = an f ads 
Y 
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Parallel 


A surface of revolution. 


PROBLEM 11-21. Instead of revolving a curve, we can revolve a two-dimensional 
region around an axis. Suppose that R is such a contented subset of the « — z plane and 
again suppose that for all (x,z) € R, x > 0. Then revolve R around the z-axis through 
360°, obtaining the solid of revolution D: 


D = {(2cos0, geme, 2) | (@,.2)e #0 <0 < 27. 


Prove that 


vol3(D) = 2m faded 
R 


2. Pappus’ theorems 


Pappus of Alexandria lived about 290-350, and has been described as the last of the great 
Greek geometers. Here is one of his theorems, which is known as the basis of modern projective 
geometry: 


26 Chapter 11 


Pappus’s theorem states: 
If the vertices of a hexagon lie 


alternately on two lines, then the 


meets of opposite sides are 


collinear. 





Just for fun, you might try proving this for yourself. 
The theorems of Pappus we want to discuss are about areas and volumes of figures of 
revolution in R*. (We shall of course prove them with calculus). 














First consider the case of a surface of revolution, using the notation introduced above. 
Then 


area(M) = an f ds. 
se 
But also the centroid of y has its x-coordinate defined by 


1, ads 
ie ds 


J, ads 
L +] 








where L = the length of the curve y. Thus we obtain Pappus’ theorem, 
area(M) = 27Z@L. 


Notice that 27% is the distance the centroid of 7 travels in one rotation. Thus Pappus’ theorem 
can be stated in words as 





The area of a surface of revolution equals the product of the length of the 
initial curve and the distance traveled around the axis by its centroid. 
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A great example is provided by a torus of revolution. Using the notation of Section 5G, where 
the center of the initial circle of radius a was located at distance b > a from the axis, the 
resulting area is given by 27a -2rb = 4r7ab. 

Sometimes the Pappus theorem can be used “backwards.” For example, suppose we revolve 
a semicircle around the z-axis as shown: 





The resulting surface is a sphere of radius a. Since we already know the area of the sphere, 
we obtain from Pappus, 
Ana* = na-2nk, 


and we have found the centroid of a semicircle: 


PROBLEM 11-22. Prove another Pappus theorem for a solid of revolution, and 
express the result in words as we did above. 














PROBLEM 11-23. Use Pappus to find the centroid of a semiball in R®. 


3. Tubes 


There is an interesting extension of Pappus’ theorems to figures in R? we might describe 
as “tubes.” Start with a nonselfintersecting C! curve 7 in R*, parametrized as y(t), a < t < b. 
Then imagine constructing a circle of radius r centered at each point y(t) of the curve, and 
orthogonal to the curve. To make sense of the orthogonality, we of course require the tangent 
vector 7/(t) # 0. Then it can be shown that if r is sufficiently small, “self-intersections” do 
not occur, and the resulting set is an actual two-dimensional manifold. We want to find its 
area. In fact, it equals the length of y times 27r. 
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In order to prove this with calculus, we need to construct an orthogonal frame at each 
point y(t). For the first vector of the frame we can use 7/(t). The other two vectors may be 
chosen rather arbitrarily. Let us call such a choice u(t), v(t), where u(t) and v(t) are unit 
vectors in R*, and the three vectors, 














y(t), w(t), v(t), 


are mutually orthogonal. We also need to know that u(t) and v(t) are of class C1. (One way to 
insure this is to assume 7¥ is of class C? and take u(t) as the unit vector in the same direction 


as : 
d y(t) 
dt | ON 





and then 


V@) 
ui) = x u(t).) 
Ily(é) | 
Now notice that differentiating the equation u eu = 1 produces u’ eu = 0. Thus we may 
express u’ as a linear combination of y’, u, and v, and the result looks like this: 





ul(t) = a(t)7(t) + e(t)u(t). 


Likewise, 
u'(t) = b(t)7(t) + d@)u(d). 
Notice also that uev = 0 implies u’ev+uev’ = 0, so that c(t) + d(t) = 0. Thus we have the 
formulas 
uo = ay’ +cv, 
vo = by — cu. 
Now we present a parametrization of our tube: 


F(t,0) = y(t) +rcosdu(t) + rsin 6v(t). 


We calculate as follows: 


OF 
aa 7 +rcos6(ay’ + cv) +r sin 6(by' — cu) 

= (1+ arcosé + brsin6)y — crsin 6u + cr cos 0v, 
OF 


— = —rsindu+rcos 6v. 


oo 
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Thus we have 


lal-7 
a0 || 
_ OF OF 
ae e 30 = cr’ sin?6+ cr? cos?@ = cr’. 
Thus the Gram determinant is therefore 
OF |2 
T(t,0)? = | OE r? — (cr’)? 








= {(1+arcosé + brsin 0)” |ly/|? + Cer Ser 
= (1+arcos@ + brsin 6)? ||y7'|/?r?. 


How nice! If r is sufficiently small that 1+ arcos@-+ brsin@ is always positive, then 
J(t,0) = (1+ arcosé + brsin@) ||9’||r. 
We conclude that 


area(M) = rf [a ateyrcose + o(e)rsiné) \|7’(t) ||dt 


b 
=r f an iv @lat 
= 2nr- length of y. 


PROBLEM 11-24. Perform the analogous computation for the 3-dimensional volume 
of a “solid” tube in R®, using the parametrization given above written now as 














F(t,6,r) = y(t)+rcoséu(t) + rsin dv(é), 
where the new variable r lies in the interval 0 < r < ro. 


4. Ellipsoids 


You probably know that the arc length of an ellipse cannot be expressed in elementary 
terms, unless of course the ellipse is merely a circle. In fact, it is for this very reason that a 
class of special functions has been defined, the elliptic integrals. Here is one example: 


DEFINITION. The complete elliptic integral of the second kind is the function EF given by 


m/2 
E(k) = | V1—k?sin*¢ dt. 
0 
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For our purposes, we restrict attention to 0 < k < 1. Of course we have the elementary values 


PROBLEM 11-25. Find one of the few known “elementary” values of EF: 


(a) ~ Vi (ath + rd) 





(see Problem 11-5). 


PROBLEM 11-26. Prove that the length of an ellipse is equal to 


4-(semimajor axis) - E(eccentricity). 











We should expect that calculating the area of an ellipsoid in R® would lead to even more 
complicated integrals. Indeed it does. However, there is quite a surprise: the case of ellipsoids 
of revolution is completely elementary. The results are given in 





PROBLEM 11-27. Consider the ellipsoid of revolution in R®, 


2 2 2 
£ Yy ae 
mpirpatp tl 


Oblate spheroid: Show that if a > b, then the area is equal to 
Qrab? loo { 2 + Ja? — b? 
Va? — b? : b , 


Prolate spheroid: Find the area in case a < b (notice that the above formula doesn’t make 
sense in this case). 


Qra? + 

















Although the length of an ellipse in R? and the area of an ellipsoid in R? are not elementary 
integrals, there happen to be associated integrals that are. Here is an example: 
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PROBLEM 11-28. Let M CR” be the ellipsoid described implicitly as 


Through any x € M construct the hyperplane 
which is tangent to M at x. Let D(x) = the 
distance from this hyperplane to the origin. 


a. Prove that 





b. Calculate explicitly 


| Ddvol,_\ : 
M 





